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ON THE X-RAY TRANSFORM OF PLANAR SYMMETRIC 

2-TENSORS 

KAMRAN SADIQ, OTMAR SCHERZER, AND ALEXANDRU TAMASAN 


Abstract. In this paper we study the attenuated X-ray transform of 
2-tensors supported in strictly convex bounded subsets in the Euclidean 
plane. We characterize its range and reconstruct all possible 2-tensors 
yielding identical X-ray data. The characterization is in terms of a 
Hilbert-transform associated with A-analytic maps in the sense of Bukhgeim. 


1. Introduction 

This paper concerns the range characterization of the attenuated X-ray 
transform of symmetric 2-tensors in the plane. Range characterization of 
the non-attenuated X-ray transform of functions (0-tensors) in the Euclidean 
space has been long known HOllTIKTOl, whereas in the case of a constant 
attenuation some range conditions can be inferred from [fTTl [H 13. For a 
varying attenuation the two dimensional case has been particularly interest¬ 
ing with inversion formulas requiring new analytical tools: the theory of A- 
analytic maps originally employed in Q, and ideas from inverse scattering 
in [[23 . Constraints on the range for the two dimensional X-ray transform 
of functions were given in [l25l |4]|. and a range characterization based on 
Bukhgeim’s theory of A-analytic maps was given in [l^ . 

Inversion of the X-ray transform of higher order tensors has been formu¬ 
lated directly in the setting of Riemmanian manifolds with boundary |j33- 
The case of 2-tensors appears in the linearization of the boundary rigidity 
problem. It is easy to see that injectivity can hold only in some restricted 
class: e.g., the class of solenoidal tensors. For two dimensional simple 
manifolds with boundary, injectivity with in the solenoidal tensor fields has 
been establish fairly recent: in the non-attenuated case for 0- and 1-tensors 
we mention the breakthrough result in Il29l . and in the attenuated case in 
031; see also lfT3 for a more general weighted transform. Inversion for the 
attenuated X-ray transform for solenoidal tensors of rank two and higher 
can be found in [|23 . with a range characterization in If23 . In the Euclidean 
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case we mention an earlier inversion of the attenuated X-ray transform of 
solenoidal tensors in [[T^ : however this work does not address range ehar- 
aeterization. 

Different from the reeent eharaeterization in terms of the scattering rela¬ 
tion in [f28l . in this paper the range eonditions are in terms of the Hilbert- 
transform for A-analytic maps introdueed in Our eharaeterization 

ean be understood as an explieit deseription of the seattering relation in 
[12^ IZTl [2^ partieularized to the Euelidean setting. In the suffieieney part 
we reeonstruct all possible 2-tensors yielding identieal X-ray data; see (l30l) 
for the non-attenuated ease and (1821) for the attenuated ease. 

For a real symmetrie 2-tensor F G 

few)- 

and a real valued funetion a G the a-attenuated X-ray transform 

ofF is defined by 

{F{xtO) 9,9) exp a{xs9)ds 

where 6* is a direction in the unit sphere S\ and (•, •) is the sealar product in 
For the non attenuated ease a = 0 we use the notation XF. 

In this paper, we eonsider F be defined on a strietly eonvex bounded set 
O C with vanishing traee at the boundary F; further regularity and the 
order of vanishing will be speeified in the theorems. In the attenuated case 
we assume a > 0 in O. 

For any (x, 6*) G O x let t{x, 9) be length of the ehord in the direction 
of 9 passing through x. Fet also consider the ineoming (—), respeetively 
outgoing (-f) submanifolds of the unit bundle restrieted to the boundary 

(3) r± ;= {(x,0) G r X : ±0-n(x) > 0}, 
and the variety 

(4) To := {(x,0) G r X : 0-ri(x) = 0}, 


(2) XaF{x,9):= 



where n{x) denotes outer normal. 

The a-attenuated X-ray transform of F is realized as a funetion on by 
/■o 

(5) X„F(x, 9) = (F(x + t9) 9,9) ^ 

J—t{x,8) 

We approaeh the range eharaeterization through its eonneetion with the 
transport model as follows: The boundary value problem 

(6) 9 ■ Vm(x, 9) + a{x)u{x, 9) = {F{x)9,9) {x, 6^) G x S^, 

(7) u\r. = 0 
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Figure 1 . Definition of /± 

has a unique solution in x and 

(8) u\r^{x,d)=X,Y{x,d), {x,9)er+. 

The X-ray transform of 2-tensors oeeurs in the linearization of the bound¬ 
ary rigidity problem ll^ : For e > 0 small, let 

g^{x) := I -f eF(x) -f o(e), x e 

be a family of metries perturbations from the Euelidean, where I is the 
identity matrix and F is as in ([T]). For an arbitrary pair of boundary points 
x,y e r let de{x, y) denote their distanee in the metrie The boundary 
rigidity problem asks for the reeovery of the metrie g" from knowledge of 
de{x, y) for all x,y e F. In the linearized ease one seeks to reeover F(x) 
from ^ |^_Q d‘^{x, y). Taking into aeeount the length minimizing property of 
geodesie one ean show that 

Id 

- - d^^{x,y)= {F{x + te)e,e)dt = xF{x,e), 

\x y\ da ^_g J—^x—y\ 

where 

\x-y\ 


2. Preliminaries 

In this seetion we briefly introduee the properties of Bukhgeim’s A- 
analytie maps [[71 needed later. 

For z = Xi + ix 2 , we eonsider the Cauehy-Riemann operators 

d = {dx, + idx^)/2, d = {dx, - idx^) 


(9) 
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Let /oo(, ^i) be the space of bounded (, respectively summable) sequences, 
£ : Zoo —^ ^oo be the left shift 


Definition 2.1. A sequence valued map 

z HA u(z) := {u_i{z),u_ 2 {z),u_^{z),...) 
is called £-analytic, ifw G C{Vt] loo) H loo) and 


9u(z) + Cd\i{z) = 0, z e Q. 


( 10 ) 


For 0 < a < 1 and Zc = 1, 2, we recall the Banach spaces in ff30l : 


OO 


(11) Z^^r) := < u = (m_i,'U_ 2 , ...) : snp'^f\u-j{C)\ < oo 



By replacing F with O and Zi with loo in (fT^ we similarly define C"(f2; Zi), 
respectively, (^“(O; loo)- 

At the heart of the theory of A-analytic maps lies a Cauchy-like integral 
formula introduced by Bukhgeim in [|3. The explicit variant (fT3l) appeared 
first in Finch [HI . The formula below is restated in terms of £-analytic maps 
as in on . 

Theorem 2.1. OTl Theorem 2.11 For some g = {g_i, g_ 2 , 9 - 3 , ■■■) £ 

C'^{F ; Zi) define the Bukhgeim-Cauchy operator B acting on g, 


n3 zv^ ((i3g)_i(z), {Bg)_2{z), {Bg)_3{z ),...), 


by 



j=0 


(13) 



Then Bg G Zoo) H C(f2; loo) and it is also C-analytic. 


For our purposes further regularity in Bg will be required. Such smooth¬ 
ness is obtained by increasing the assumptions on the rate of decay of the 
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terms in g as explicit below. For 0 < a < 1, let us recall the Banach space 
Ya in CT : 


(14) Y^={ge e(r) : sup < oo 






Proposition 2.1. |l3ll Proposition 2.1] Ifg e « > 1/2, then 

(15) Bg e h) n ^"(n; h) n loo). 

The Hilbert transform associated with boundary of ^-analytic maps is 
defined below. 


Definition 2.2. Forg = {g_x,g_ 2 ,g-z,...) G l)^{r)r\C°‘{r-,li), we define 
the Hilbert transform Fig componentwise for n > 1 by 


(Hg)_„K) = 


TT 


(16) 


+ 


9-niC) 

c-e 

dC 


dc 


dC 


^ Jr [C - ^ (- U 


} oo 

J2d-n-ji0 


c-e 

c-e 


, eer. 


The following result justifies the name of the transform FL. For its proof 
we refer to If^ Theorem 3.2]. 


Theorem 2.2. For 0 < a < 1, let g e lU{T) n C'"(r;/i). For g to be 
boundary value of an C-analytic function it is necessary and sufficient that 

(17) (J + z7()g = 0 , 

where Ft is as in CH). 


3. The non-attenuated case 


In this section we assume a = 0. We establish necessary and sufficient 
conditions for a sufficiently smooth function on T x SMo be the X-ray 
data of some sufficiently smooth real valued symmetric 2-tensor F. For 
6 = (cos 99 , sin p) G S^, a calculation shows that 

(18) (F(x)0, 9) = fo{x) + , 


where 

(19) /o(a;) =-^-, and f 2 {x) 

The transport equation in ® becomes 


fll{x)-f22{x) , .fl2{x) 

4 2 


(20) 9 ■ Vu{x, 9) = fo{x) + f 2 {x)e^'^ + x e n. 
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For z = xi 1x2 G we consider the Fourier expansions of u{z,-) in 
the angular variable 6 = (cos ip, sin ip)\ 

OD 

— 00 

Since u is real valued its Fourier modes occur in conjugates, 

U-n{z) = Un{z), n, > 0, z E VL. 

With the Cauchy-Riemann operators defined in dH) the advection operator 
becomes 

e-V = e-^^d + e^^d. 


Provided appropriate convergence of the series (given by smoothness in the 
angular variable) we see that if u solves (l20l) then its Fourier modes solve 
the system 

(21) dui{z) + du-i{z) = fo{z), 

(22) du_i{z) + du_ 3 {z) = f 2 {z), 

(23) du 2 n{z) + du 2 n- 2 {z) = 0, n < 0, 

(24) dU 2 n-l{z) + du2n-3{z) =0, n< -1, 

The range characterization is given in terms of the trace 

(zjj g u\rxs^ — 


0 , {x, 6) e r_U Fq. 

More precisely, in terms of its Fourier modes in the angular variables: 

CX> 

(26) s(C.«) = y^Sn(C)e‘"*’, (er. 

— OD 

Since the trace g is also real valued, its Fourier modes will satisfy 


(27) g-niC) = gniO, n>0, Cer. 

From the negative even modes, we built the sequence 

(28) g^'’^-:={go,g-2,g-,,...). 

From the negative odd modes starting from mode —3, we built the sequence 


(29) 


g° ■—{g-3, g-5, g-7, ■■■)■ 


Next we characterize the data g in terms of the Hilbert Transform "H in 
([HI). We will construct simultaneously the right hand side of the transport 
equation (l20l) and the solution u whose trace matches the boundary data 
g. Construction of u is via its Fourier modes. We first construct the neg¬ 
ative modes and then the positive modes are constructed by conjugation. 
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Except from negative one mode m_i all non-positive modes are defined by 
Bukhgeim-Cauehy integral formula in (fT3l) using boundary data. Other then 
having the traee g_i on the boundary m_i is unconstrained. It is chosen ar¬ 
bitrarily from the elass of funetions 

(30) := jv' e C) : i/;|r= ^7-i} • 


Theorem 3.1 (Range eharacterization in the non-attenuated case). Let a > 

1 / 2 . 


(i) Let F G 


^1,0 

^0 


(0;M' 


2x2^ 


For g : = 


consider the corresponding sequences as in 
Then g°'^'^ e (^^(r) fl C“(r; /i) satisfy 


XF{x,9), {x,9)er+, 

0 , {x, 9 ) E r_ u Fq, 

as in 


and 


(31) [I + = 0, 

(32) [/ + = 0, 

where the operator Ft is the Hilbert transform in (fT^ . 

(ii) Let g G C" (F; C^’“(S^)) fl C'(F; C^’"(S^)) be real valued with 
5'|r_uro= 0. If the corresponding sequence g^^^en godd ^ satisfies (l3TI) 
and (l32l) . then there exists a real valued symmetric 2-tensor F G F(O; 
such that g\r^= XF. Moreover for each E in (l30l) . there is a unique 
real valued symmetric 2-tensor such that g\r+= XF.^. 


Proof (i) Necessity 

Let F G Fd’“(f2; Since F is eompactly supported inside Q, for 

any point at the boundary there is a eone of lines whieh do not meet the 
support. Thus = 0 in the neighborhood of the variety Fq which yields g E 
F^’"(F X S^). Moreover, g is the traee on F x of a solution u E F^’“(f2 x 
S^) of the transport equation (l20l) . By [1^ Proposition 4.1] g 

e(r)nF“(r;/i). 

If u solves (l20l) then its Fourier modes satisfy (l2T]) . (l22l) . (l23l) and (l24l) . 
Sinee the negative even Fourier modes U 2 n of u satisfies the system (|2^ for 
n < 0, then 


^ ^ u®^®"(2;) := {uo{z),U-2{z),U-fiz),U-Q{z), ■■■) 

is F-analytie in O and the neeessity part in Theorem |2]^ yields (ISTl) . 

The equation (l24l) for negative odd Fourier modes U 2 n-i starting from 
mode —3 yield that the sequence valued map 

z H- u^'^{z) := {u_3{z),u_5{z),u_7{z), ■■■) 

is F-analytic in fl and the necessity part in Theorem [2]^ yields (l3^ . 

(ii) Sufficiency 
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To prove the suffieieney we will eonstruet a real valued symmetrie 2- 
tensor F in 12 and a real valued funetion u G C^(12 x S^) fl x S^) sueh 
that M|rxsi= 9 u solves (l20l) in 12. The eonstruetion of sueh u is in 
terms of its Fourier modes in the angular variable and it is done in several 
steps. 

Step 1: The construction of negative even modes U 2 n for n < 0. 

Let (7 G C" (r';C^’“(S^))nC'(r';C'^’"(S^))be real valued with 5(|r_uro = 

0 . Let the eorresponding sequenees satisfying (ISTl) and satisfying 

(l32l) . By [1^ Proposition 4.1(ii)] goda ^ Bukhgeim- 

Cauehy Integral formula (fTSl) to eonstruet the negative even Fourier modes: 

(33) {uo{z),U-2{z),U-i{z),U-e{z ),...) := z G 12. 

By Theorem 12.11 the sequenee valued map 

X H- {uo{z),U-2{z),U-4{z), ...), 

is £-analytie in 12, thus the equations 

(34) du-2k + du-2k-2 = 0, 

are satisfied for all A; > 0. Moreover, the hypothesis (l3TI) and the suffieieney 
part of Theorem I2.2l vields that they extend eontinuously to F and 

(35) u_2k\r = 9-2k, k>Q. 

Step 2: The construction of positive even modes U 2 n for n > 1. 

All of the positive even Fourier modes are eonstrueted by eonjugation: 

(36) U2k ■= k>l. 

By eonjugating (l34l) we note that the positive even Fourier modes also sat¬ 
isfy 

(37) du2k+2 + du2k = 0, k>0. 

Moreover, they extend eontinuously to F and 

(38) U2k\r = w_2fc|r = g-2k = 92k, k>l. 

Thus, as a summary, we have shown that 

(39) du2k + du2k-2 = 0, Vfc G Z, 

(40) U2k\r= 92k, VA: G Z. 

Step 3: The construction of modes m_i and ui. 

Let G 4/g as in (l30l) . We define 

(41) U-i := Ip, and ui := xp. 

Sinee g is real valued, we have 

(42) ui\r=WF = 9i- 
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Step 4: The construction of negative odd modes U 2 n-i for n < — 1. 

Use the Bukhgeim-Cauchy Integral formula ([13]) to construct the other 
odd negative Fourier modes: 

(43) (u_3(^),m_5(^), •••);= zeVL. 

By Theorem l2.1[ the sequence valued map 

2 ; H- {U-z{z),U-^{z),U-7{z), 
is ^-analytic in U, thus the equations 

(44) du2k-i + du2k-3 = 0, 

are satisfied for all A: < —1. Moreover, the hypothesis (l3^ and the suffi¬ 
ciency part of Theorem 12.21 yields that they extend continuously to F and 

(45) U2k-i\r = g2k-i, VA: < —1. 

Step 5: The construction of positive odd modes U 2 n+i for n > 1. 

All of the positive odd Fourier modes are constructed by conjugation: 

(46) U2k+3 ■= M-( 2 fc-r 3 ), A: > 0. 

By conjugating (l44l) we note that the positive odd Fourier modes also satisfy 

(47) du2k+3 + du2k+i =0, VA: > 1. 

Moreover, they extend continuously to F and 

(48) U2k+3\r = U_{2k+3)\r = g-{2k+3) = g2k+3, k > 0. 

Step 6: The construction of the tensor field whose X-ray data is 

g- 

We define the 2-tensor field 

/o -f 2 Me /2 
2 Im /2 

where 

(50) /o = 2Me(c)-0), and /2 = 

In order to show g\r+= XF^ with as in (|49|) . we define the real valued 
function u via its Fourier modes 

(51) u{z, 9) := uq{z) -f ip{z)e~'''^ + 'tp{z)F'^ 

oo oo 

n=2 n=2 

and check that it has the trace g on F and satisfies the transport equation 


(49) 


:= 


2Im/2 
/o — 2 Me /2 


Since c/ G U" (F; C^’^{S^))nC{F-, ^^’“(Si)), we use (SOl Corollary 4.1] 
and Proposition 4.1 (iii)] to conclude that u defined in (|5TI) belongs to 
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X S^) n C°‘{VL X S^). In particular m(-,6*) for 6 = (cos(^, sin(^) 
extends to the boundary and its traee satisfies 


u{-,9)\r= iuQ + ipe 


Un^ 


in^p 


n=2 


n=2 


— Wo|r+'0|re + "^IrM-n|re Mn|re' 


imp 


n=2 


n=2 


— fi'o + 9-1^ S' 


— 


—inp 


n=2 


+ ^9r 

n=2 


pinp 


where in the third equality above we used (l40l) . (I45l).(l48l). (l42l) and definition 
of V' e in (1^ . 

Sinee m G x S^) fl C"(f2 x S^), the following ealeulation is also 

justified: 


Q ■ Vu = e + e^'^duo + +e + d%jj + d%jj + 

OO OO 

n=2 n=2 

OO OO 

n=2 n=2 

Rearranging the modes in the above equation yields 

6 ■ Vm = (d'll) + du-s) + e^'^'^{d'ip + du^) + dip + d'lfj 
+ e~'‘‘^{duo + du- 2 ) + e*‘^(9Mo + du2) 

OO OO 

n=l n=l 

Using (1^ . (HH), and (l47l) simplifies the above equation 

S' ■ Vm = + £^*‘^(5-0 + cIms) + clV' + 

Now using (I50l) . we eonelude (l20l) . 

e-Vu = e-2*^/2 + + /o = (Fv,0,0). 


□ 

As the source is supported inside, there are no ineoming fluxes: hence the 
traee of a solution u of (l20l) on r_ is zero. We give next a range condition 
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only in terms of g on r+, where g := M|rxsi- More preeisely, let u be the 
solution of the boundary value problem 


9 ■Vu{x,e) = {F{x)9,9), xen, 

^ ^ Hz, 9) =-^g\r+{z,-9), {z,9) e r_. 

Then one ean see that 

(53) u\r+= -^g\r+, 

and therefore wlrxsi is an odd funetion of 9. This shows that we ean work 
with the following odd extension: 


(54) 


9{z,9) 


9{z,9) -g{z,-9) 
2 


{z,9)e{rxS^)\ro, 


and ^ = 0 on Fq. Note that g is the traee of m on T x 
The range eharaeterization ean be given now in terms of the odd Fourier 
modes of g, namely in terms of 


(55) 


g {9-3, 9-5, 9-7, ■■■)■ 


Corollary 3.1. Let a >1/2. 

(i) Let F G Cg’“(f2; u be the solution of (l52l) and g as in (l55l) . 

Then g G /^^(T) fl (7“(r; li) and 

(56) [I + i'Hjg = 0, 

where the operator Ft is the Hilbert transform in (fT^ . 

(ii) Let g G C" (F; C^’“(S^)) fl (^(F; C^’"(S^)) be real valued with 
5'|r_uro= 0. Let g be its odd extension as in (l54l) and the corresponding 
g as in (l55l) . If g satisfies (l5^ . then there exists a real valued symmet¬ 
ric 2-tensor F G F(f2;M^^^), such that g\r^= XF. Moreover for each 
i/ G 4/g in (l30l) . there is a unique real valued symmetric 2-tensor F^ such 
that g\r+= XF^. 


4. The attenuated case 

In this seetion we assume an attenuation a G C'^’“(f2), a > 1/2 with 

nun a > 0. 
n 

We establish neeessary and suffieient eonditions for a suffieiently smooth 
funetion on F x SMo be the attenuated X-ray data, with attenuation a, 
of some suffieiently smooth real symmetrie 2-tensor, i.e. g is the traee on 
F X of some solution u of 

(57) 9-Vuix^b)+ a{x)u{x,9) = {F{x)9,9), {x,9)erxS\ 
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Different from 1-tensor ease in [l3T1l (where there is uniqueness), in the 
2-tensor ease there is non-uniqueness: see the elass of funetion in (|8^ . 

As in [f30l we start by the reduetion to the non-attenuated ease via the 
speeial integrating faetor e~^, where h is explieitly defined in terms of a by 

(58) h{z, 6) := Da{z, 0) - ^ (/ - iH) Ra{z ■ 6^, 6), 


where 6-^ is orthogonal to 6, Da{z, 0)= a{z -f t6)dt is the divergenee 

Jo 

POO 

beam transform of the attenuation a, Ra{s, 0)= a (sd'^ + t6) dt is the 


Radon transform of the attenuation a, and the elassieal Hilbert transform 

1 r h{t) 


Hh{s) = - 


-dt is taken in the first variable and evaluated at s = 


TT 


t 


z- 6-^. The funetion h was first eonsidered in the work of Natterer [|2T]I : see 
also [O, and [jH for elegant arguments that show how h extends from 
inside the disk as an analytie map. 

The lemma 1431 and lemma l43] below were proven in [|3TI for a vanishing 


at the boundary, a G Co’“(f2), a > 1/2. We explain here why the vanishing 
assumption is not neeessary: we extend a in a neighbourhood (2 of H with 
eompaet support, d G Cq’"(H). We apply the results (STl Lemma 4.1 and 
Lemma 4.2] for the extension d and use it on H. 


Lemma 4.1. [|3T1 Lemma 4.1] Assume a G C^’’"(D), p = 1,2, a > 1/2, 
and h defined in (l58]) . Then h G C^’“(H x S^) and the following hold 
(i) h satisfies 


(59) 9 ■ Vh{z, 9) = -a{z), {z, 9) e fl x S\ 


(ii) h has vanishing negative Fourier modes yielding the expansions 
(60) 

oo oo 

:= ^ , (x, 0) G H x S\ 

fc=0 k=0 

with 

(Hi) 

(61) z HA {afiz), a 2 (z), 03(2:), ...,) G CP’"(f2; /i) n ^(fl; If, 

(62) z ^ {I3fz),l32{z),l3fz), ...,) G h) n C(n- If. 
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(iv) For any z 

(63) dPo{z) = 0, 

(64) dl3i{z) = -a{z)l3o{z), 

(65) dl3k+2{z) + dl3k{z) + a{z)l3k+i{z) = 0, fc > 0. 

(v) For any z eFL 

( 66 ) dao{z) = 0 , 

(67) dai{z) = a{z)ao{z), 

( 68 ) dak+ 2 {z) + dak{z) + a{z)ak+i{z) =0, k > 0. 

(vi) The Fourier modes ak, l3k,k > 0 satisfy 

k 

(69) 1, ^ ^ Oim/^k—m 0, fc ^ 1. 

m=0 

From (l59l) it is easy to see that u solves (l57l) if and only if u := e~^u 
solves 

(70) e ■ Vv{z, 6) = {F{z)e, 

If u{z, 9) = J2'^=-oo Un{z)e^'^‘^ solves (l57l) . then its Fourier modes satisfy 

(71) dui{z) + du_i{z) + a{z)uQ{z) = fo{z), 

(72) duo{z) + du_ 2 iz) + a{z)u_i{z) = 0, 

(73) du-i{z) + du- 3 {z) + a{z)u- 2 {z) = f 2 {z), 

(74) dun{z) + dun- 2 iz) + a{z)un-iiz) =0, n < -2, 

where /o, /2 as defined in (fT9l) . 

Also, if V := e~^u = X]^-oo solves (TTOl) . then its Fourier 

modes satisfy 

dvi{z) + dv_i{z) = ao{z)fo{z) + a 2 iz)f 2 iz), 
dvo(z) + dv- 2 {z) = ai{z)f 2 {z), 
dv-i{z) + dv- 3 {z) = ao{z)f 2 {z), 

(75) dvniz) + dvn- 2 iz) = 0, n< -2, 

where a^, Oi and 02 are the Fourier modes in (l60l) . and /o, /2 as defined in 

dH. 

The following result shows that the equivalenee between (1741) and (1751) is 
intrinsic to negative Fourier modes only. 
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Lemma 4.2. [[311 Lemma 4.2] Assume a G a > 1/2. 

(i) Letv = {v- 2 , V- 3 ,...) G /i) satisfy (1751) . and u = (m_ 2 , M-s, •••) 

be defined componentwise by the convolution 

OO 

(76) Un := n < -2, 

i=o 

where fij ’s are the Fourier modes in (l60l) . Then u solves (1741) m 0. 

(ii) Conversely, let u = (m_ 2 , M-s, ■••) G If) satisfy (1741) . and v = 

(u_ 2 , V- 3 ,...) defined componentwise by the convolution 

OO 

(77) Vn ■.= '^ajUn-j, n<-2, 

j=0 

where aj’s are the Fourier modes in (l60l) . Then v solves (1751) in (2. 

The operators 0, 9 in dH) can be rewritten in terms of the derivative in 
tangential direction and derivative in normal direction dn, 

dn = COS + sin , 
dr = — sin Tjdx^ + cos , 

where p is the angle made by the normal to the boundary with xi direction 
(Since the boundary F is known, p is a known function on the boundary). 
In these coordinates 

-ir] _ ir) 

(78) d = ^^{dn-idr), d=—{dn + idr). 

Next we characterize the attenuated X-xay data g in terms of its Fourier 
modes go, g-i and the negative index modes 7 _ 2 , 7 - 3 , 7 - 4 ... of 

OO 

(79) 7fc(C)e*''^ C ^ T. 

k=—oo 

To simplify the statement, let 

(80) gfe := (7_2,7-3,7-4..-), 

and from the negative even, respectively, negative odd Fourier modes, we 
built the sequences 

(81) gr" = (7-2,7-4,...), and = ( 7 - 3 , 7 - 5 , •••)• 

Note that 7 _i is not included in the definition. As before we construct 
simultaneously the right hand side of the transport equation (l57l) together 
with the solution u. Construction of u is via its Fourier modes. We first 
construct the negative modes and then the positive modes are constructed 
by conjugation. Apart from zeroth mode uq and negative one mode m_i, 
all Fourier modes are constructed uniquely from the data g/'^'^. The 
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mode Mo will be ehosen arbitrarily from the elass with preseribed traee 
and gradient on the boundary F defined as 


(82) V>eC'2(f];M):V>|r=^o, 


dn'4’\r= —2Mee 


+a|rc/-ij 




where B be the Bukhgeim-Cauehy operator in (fT3l) . (3/s are the Fourier 
modes in (1^ and gh in (l80l) . The mode m_i is define in terms of Uq, see 


Reeall the Hilbert transform FL in (fT^ . 


Theorem 4.1 (Range eharaeterization in the attenuated ease). Let a e (H), 

a > 1/2 with nfin a > 0. 
n 


(i) Let F G For g : = 

consider the corresponding sequences 
l]^{T)FC-{T-h) satisfy 


fX,F(a:,0), (a;,0)Gr+, 

\ 0, {x, 0) G F_ U Fq, 

^even^ ^odd ^ ^even ^ ^odd 


G 


(83) [/ + *-H]gr” = 0, [/ + *-H]gr = 0, and 


(84) dr go = -2Ime | d'^ f3j{Bgh)-2-j 

j=0 


+ o.\r 9-1 , 


where Ft is the Hilbert transform in (fT^ . B is the Bukhgeim-Cauehy oper¬ 
ator in (fT3l) . ’s are the Fourier modes in (l60l) and gh in (f80l) . 

(ii) Let g G C“ (F; C'^’“(S^)) fl F(F; F^’"(S^)) be real valued with 
fi'|r_uro= 0- If the corresponding sequences , gf^'^ G satisfying 
([83]) and (l84l) then there exists a symmetric 2-tensor F G F(f2; such 

that g\r+= XaF. Moreover for each i/ G in (l82l) . there is a unique real 
valued symmetric 2-tensor such that g\r_^.= XaF^. 


Proof (i) Necessity 

Let F G Fd’“(f2; Sinee F is eompaetly supported inside Ft, for 

any point at the boundary there is a eone of lines whieh do not meet the 
support. Thus = 0 in the neighborhood of the variety Fq whieh yields 
g G (7^’"(F X S^). Moreover, g is the traee on F x of a solution u G 
Fi’"(n X S^). By lllOl Proposition 4.1] gf'^ G llf{T) n F"(r; h). 

Let V := e~^u = , then the negative Fourier modes 

of V satisfy (fTSl) . In partieular its negative odd subsequenee (m_ 3 , m_ 5 , ...) 
and negative even subsequenee (m_ 2 , M- 4 , ...) are £-analytie with traees gf^ 
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respectively The necessity part of Theorem 12.21 yields (l8^ : 

[/+*?^]gr = 0, [/+^?^]gr" = o. 

If u solves (l57l) . then its Fourier modes satisfy (ItT]) . (17^ . (1731) . and (1741) . 
The negative Fourier modes of u and v are related by 

OO 

(85) Un = '^(3jVn-j, n < 0, 

j=0 

where /9/s are the Fourier modes in (l60l) . The restriction of (1721) to the 
boundary yields 

duo\r = —du-2\r—{au-i)\r- 

Expressing d in the above equation in terms of dr and dn as in (1781) yields 

— + idr)uo\r = -du_ 2 \r-a\rg-i- 

Simplifying the above expression and using drUo\r= dr go, yields 
dnUo\r+idrgo = {du-2\r+a\rg-i) ■ 

The imaginary part of the above equation yields (IMl) . This proves part (i) 
of the theorem. 

(ii) Sufficiency 

To prove the sufficiency we will construct a real valued symmetric 2- 
tensor F in and a real valued function m G x S^) D (7(f2 x S^) such 
that M|rxsi= g u solves (1571) in fl. The construction of such u is in 
terms of its Fourier modes in the angular variable and it is done in several 
steps. 

Step 1: The construction of negative modes Un for n < —2. 

Let^f G (F; (7^’“(S^))n(7(r'; (7^’“(S^)) be real valued with 5(|r_uro = 

0. Let the corresponding sequences in (IMl) satisfying (l83l) and 

dHU). By |l30l Proposition 4.1(ii)] and [l30l Proposition 5.2(iii)] g^®”, G 
Ya. Use the Bukhgeim-Cauchy Integral formula (fT3l) to define the ^-analytic 
maps 

(86) v--(x) = {v_2{z),v_,{z ),...) ;= i3gr"(^), e (2, 

(87) = {v.o{z),v_,{z), ...) := i3gf"(z), x G U. 

By intertwining let also define 

v(2;) := (r;_2(2:),v_3(2:),...), z E il. 

By Proposition 12. II 

(88) V G h) n ^“(n; h) n /oo). 
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Moreover, sinee satisfy the hypothesis (l8^ . by Theorem 12.21 we 

have 


^odd 


even\ _ even j odd\ _ ot 

V \r— Sh V \r— Sh 

In partieular 

OO 

(89) Vn\r='^iak\r) gn-k, n <-2. 

k=0 

For eaeh n < —2, we use the eonvolution formula below to eonstruet 

OO 

(90) Un ■ 

j=0 

Sinee a G C^’“(f2), by (1^ . the sequenee z i-)- {l3o{z),/3i{z),/32{z),...) is 
in C^’"(f2; li) fl C"(f2; h). Sinee eonvolution preserves h, the map is in 

(91) 2 ; I-)- {u_2{z),u_s{z),...) e c^’“((2; h) fl h). 

Moreover, since v G C^(f2; loo) as in ([88]), we also conclude from convolu¬ 
tion that 

(92) 2 ; H- {u_2{z),u_3{z ),...) G loo). 

The property ([93 justifies the calculation of traces Un\r for each n < —2: 

OO 

Un\r = ''^Pj\r{vn-j\r) ■ 

j=0 

Using (l89l) in the above equation gives 

OO OO 

Un\r = Efti cifclrfi'n-j-fc- 

j=0 k=0 

A change of index m = j + k, simplifies the above equation 

OO m 

^ri|_r ^ ^ ^kPm—kQn—m'j 

m=0 k=0 

00 m 

= aoPogn + EE akl3 

m—kQn—m' 

m=l k=0 

Using Lemma l4~T] (vif yields 

(93) Un\r = Qm n< -2. 

From the Lemma lL2l the constructed Un in (l90l) satisfy 

(94) dUn + dUn-2 + UMn-l = 0, tl < -2. 

Step 2: The construction of positive modes for n>2. 
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All of the positive Fourier modes are eonstrueted by eonjugation: 

(95) Un := n>2. 

Moreover using (|9^ . the traees Un\r for eaeh n > 2: 

(96) ^n|_r U—n\r 9—n Qrii n'^2. 

By eonjugating (|9^ we note that the positive Fourier modes also satisfy 

(97) dun +2 + dun + aUn+i = 0, n>2. 

Step 3: The construction of modes mq, M-i and ui. 

Let xjj G as in (f82l) and define 

(98) Mo := -0, 


and 

(99) 


M_i := 


—d'ljj — du 2 


Ml := M_i. 


By the eonstruetion mq G C^(f2; Zoo) and m_i G Zoo), and 

(100) duo + (9m_ 2 + aM_i = 0 


is satisfied. Furthermore, by eonjugating (1 1001) yields 

(101) duo + du2 + aui = 0. 


Sinee -0 G \I/“, the traee of mq satisfies 
(102) 'Uo|r= go- 


We eheek next that the traee of m_i is 


M_i r = 


-00 — c)m 2 


r 


1 

a 

1 

2a 


o^V 


-{dn + id, 


T)^g\r- 


r 


c)m_ 


2 r 


e*’'{c)„0|r+i0r0|r+2e "'^du-2\r} 


r 


(103) 


= 9-1, 


where the last equality uses (f84l) and the eondition in elass (1^ . 

Step 4: The construction of the tensor field whose attenuated X- 
ray data is g. 

We define the 2-tensor 


( 104 ) 


//o + 2 Me 02 2 Im 02 \ 

2Im02 0o-2Me02y'’ 
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where 


(105) 

(106) 


/o = -2 Me 


dif] + du-2 


+ a-^, and 


/2 = -d 


— f d%l) du-2 


+ c)m_3 + au_2- 


Note that /2 is well defined as m _2 G loo) from (|92|) . 

In order to show g\r+= with as in (11041) . we define the real 

valued funetion u via its Fourier modes 

(107) u{z, 9) := uo{z) + + ir^{z)e^‘^ 

oo oo 

n=2 n=2 

We eheek below that u is well defined, has the traee S' on F and satisfies the 
transport equation (l57l) . 

For eonvenienee eonsider the intertwining sequenee 

u( 2 ;) := {uq{z),U-i{z),U- 2 {z),U- 2 .{z), ...), z eVt. 

Sineeu G F^’“(f2;/i)nF"(f2;/i), by [SOI Proposition 4.1 (iii)] weeonelude 
that u is well defined by (11071) and as a funetion in (7^’"(f2 x S^) D (7"(f2 x 
S^). In partieular m(-, 6) for 6 = (cos ^9, sin tp) extends to the boundary and 
its traee satisfies 


u{-,9)\r= (mo + M-iC 

n=2 

oo 

= Mo|r+M-i|re“*‘^ + 

n=2 


—iriip 




Tlri.C 


m^p 


n=2 


r 




— 2 n<p 




'^nlr 


inip 


n=2 


= 9o + 9-ie ^ 9-ne 

n=2 

= 9{;0), 


—in(p 


+ E 


9n^ 


in(p 


n=2 


where is the third equality we have used (|9^ . (I96l) . (11021) . and (11031) . 
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Since u G x S^) n C"(n x S^), the following ealeulation is also 

justified: 

9 ■ Vu + au = + e^'^duo + e~‘^^'^du-i + dui + du-i + e^^'^dui 

OO OO 

n=2 n=2 

OO OO 

n=2 n=2 

OO OO 

+ auo + aM_ie“*‘^ + awie*'^ + aUnC^"''^. 

n=2 n=2 

Rearranging the modes in the above equation yields 

6 ■ Vu + au = + 9m_3 + au_ 2 ) + + 5m3 + 0 ^ 2 ) 

+ e~^'^{duo + 9'U_2 + aM_i) + e^‘^{duo + du2 + awi) 

OO 

+ dui + (9k_i + auQ + ^(c?ti„+2 + 

OO 

+ 'Y^{dU-n + dU-n-2 + 

n=2 

Using (|9^ . (l97l) . (1 1001) and (IIOII) simplifies the above equation 

9 ■ Vu + au = + 9m_3 + au_ 2 ) + e^'^‘^{dui + 5^3 + 0 ^ 2 ) 

+ (9mi + 9m_i + auQ . 

Now using (11051) and (11061) . we eonelude (l57l) 

9 ■Vu + au = e-^^^f 2 + €^^^2 + /o = (F^^, 9). 

□ 
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